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Abstract 

We study the effect of Bit continuum in the QCD sum rule analysis of the heavy 
meson doublet (0 + ,l + ) in the leading order of heavy quark effective theory. New 



[ sum rules are derived for the leading order binding energy A + i and pionic coupling 

constant q' . 

PACS number: 12.39.Hg, 11.30.Rd, 13.20.He, 12.38.Lg 

oo 
o\ , 

1 Introduction 

i ■ 

The QCD sum rules Jl|] for the masses of the excited heavy meson doublet (B' , B[) of spin 
^ . parity (0 + , l + )(+ .i) have been studied in || |3|, [|, |E|, where the indices (+, |) denote the 
parity and spin of the light component jj. Recently the 0(a s ) correction to the m(0 + ) 
sum rule has been calculated in || |7j. In ||, || the masses of the (0 + , l + )( + i), together 
with those of the doublet (1 + , |) were calculated up to the order of 0{l/rriQ) in the 
framework of heavy quark effective theory (HQET) ||. Let Apj, = Tnpj l — m&, where 
mpj l is the mass of the the doublet in the leading order of 1/tuq and P,ji are the parity 
and spin of the light component of the heavy mesons in the doublet. The results in || in 
the leading order of a s are A + i = (1.15 ± 0.10)GeV if the usual interpolating current of 
the lowest dimension is used and A + § = (0.82 ± 0.10)GeV. Within errors the results in 

0, [3], |] are consistent with these results. This would imply that + state lies 100 300 

MeV above 2 + state. The O^I/tuq) corrections calculated in || does not change much 
this mass difference. This result is inconsistent with the new experimental data where 
m(2 + ) is about 100 MeV larger than m(0 + ). It is unlikely that 0(a s ) corrections can 
account for this discrepancy. 



Recently Blok et al. made the following observation [1C]. Due to the 5- wave nature 
of the Bit intermediate state and the large coupling of the soft Goldstone particle, the 
contribution of the Bit continuum to the spectral density in the correlator of two + 
currents is unusually large. It rises faster than the quark-gluon spectral density in the low 
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energy region and exceeds it in magnitude in that region. Thus, it may violate the naive 
quark-hadron duality if we integrate the spectra over a region below the lowest pole. It 
was proposed in [TIJ that this is the reason for the abnormal large value of the residue of 
the pole JH obtained in the standard "lowest pole plus parton-like continuum model" for 
the QCD sum rules. Therefore, a better approach is to include the Bit continuum in the 
soft pion region in addition to the lowest B' pole in the sum rule. 

In this note we shall study in detail this effect in the sum rules for A + 1 and the residue 
of the B' Q pole. Besides, the decay rates of B' Q and B[ have been studied in [UJ with 



the light cone QCD sum rules (LCQSR) [0. In view of the duality violation effect due to 
the intermediate states of Bn continuum in the sum rule we shall also re-investigate the 
LCQSR for pion coupling of the (0 + , 1 + ) doublet including the contribution of Bn states. 
Section [|is a short review of previous sum rules. New QCD sum rules with intermediate 
state contribution and the numerical analyses are presented in section [| The last section 
is a short summary. 

2 Previous sum rules 

2.1 Previous mass sum rules 

The interpolating current for the doublets (0 + , 1 + ) reads 

J l,+,\ = M , (i) 

41 i = Kl"l?q , (2) 

where h v (x) is the heavy quark field in HQET, is the heavy hadron velocity, 7™ = 
7" — vv a , the indices j, +,ji in Jj t + t j t are the total angular momentum, the parity and 
the light component angular momentum respectively. Note there is a factor of 4^ in the 
definitions of the interpolating currents in our previous work. 
We consider the correlator 

II(ui) =ij rf 4 xe^(0|T{J 0i+j ,(x), Jj |+f4 (0)}|0> (3) 
and define the overlapping amplitude / + i as 

(0|J 0i+) i(0)i^) = / +) i. (4) 
Assuming the standard "pole plus parton-like continuum" model we get 

1 

II(co>) = - j h continuum . (5) 

2 A 1 1 — uj 

Here to = v ■ k, which is a factor | smaller than the u used in 0. We derive the sum rule 
for the (0+, 1+) doublet [| 

where (qq) = (qga^G^q). 
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Using the following standard values for the condensates 

(qq) = -(0.225 GeV) 3 , 
(a s GG) = 0.038 GeV 4 , 

m 2 = 0.8 GeV 2 , (7) 

and with s = (1.5 ± 0.1) GeV, T = 0.4 ~ 0.6 GeV we get from (§) 

A +j i = (1.15 ±0.10) GeV, (8) 

f +1/2 = (0.570 ± 0.08) GeV 3/2 . (9) 
For comparison here we also write down the sum rules for the (1 + , 2 + ) doublet. 

\n^ ,T - h S> 1 ^ - h < <™> ~ t^ t • < 10 > 

where the following interpolating currents for the (1 + , 2 + ) doublet are used 

4 a +i§1 = yf h^^iv?- 1 -^ p t y, (ii) 

(— i) / 2 



4T 2 = K^^^r+^r-^gr'Vtjq. m 

These currents are also a factor \/2 larger than those in || . 

2.2 Previous sum rules for pionic couplings 

Let us define the decay amplitudes of the doublet (0 + , 1 + ) in full QCD 

m? n i — m 2 R 

M(B' ^B(k)7r(q)) = Ijm^m ° # 9 , (13) 

B o 

M(B[^B*(k)n(q)) = I{e* ■ V9Jm B[ m B ^ + (k ■ v )(q ■ e*)F} , (14) 

where I = \/2, 1 for charged and neutral pion respectively. The structure F vanishes in 
the rriQ — > oo limit. 

In HQET these amplitudes have the simple form 

M(B' ^Bn) = Ig', (15) 

M(B[ -> B*n) = I e* ■ rig' (16) 

where 

gi = -g(A +h -A_ h ). (17) 

For deriving the sum rules for g' we consider the correlator 

d 4 xe- ik - x (n(q)\T{J ^ l2 (0), jj i(a;)}|0) = I G B , B (u,u;') , (18) 
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where k' = k + q, u = v ■ k, u' = v ■ k', q 2 = and Jl 1 = h v ^q. 

u, , 2 

Again after invoking the naive quark-hadron duality we have the double dispersion 
relation: 

c c' . . 

-7= ^--7= r + T + T h continuum . (19) 

4(A_i -u')(A + i -uj) A_i-a/ A +) i-w v ' 



Expressing (jig) with the pion wave functions |12] we obtain the LCQSR for g' [|TT 



= 2F.e^ : V ±A {-^(n )T 2 /i(|) + ^p(«oWo(^) + 9> )} , (20) 

where = 93MeV, //^ = — = 1.32GeV at the scale // = lGeV, ipp(u) etc are the 
light cone pion wave functions defined by \V2\ 

< 7r(g)|J(x)i7 5 w(0)|0 >= v^/V C due iuqx VP {u) . (21) 

Jo 

(p' n (u ), g'^Uo) are the first derivatives of (p n (u), gi(u) at u = u , u = T + r , T = ^pfr, 
Ti, T2 are the Borel parameters. We choose T\ = T<i- Therefore uq = \. At this point 
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(p'(uo) and g[(u ) vanish. The factor f n (x) = 1 — e~ x fr is used to subtract the 

parton-like continuum contribution with the continuum threshold Sq. 
Correcting a numerical error in || we obtain from (|20"D 

^ = 3.6 ±0.7. (22) 



3 New sum rules with Bn intermediate states 

3.1 Mass sum rules 

We have the dispersion relation for (|]) 




where p(s) is the spectral density in the limit ttiq — > 00. At the quark level, 

P 9 (s) = ^s 2 , (24) 

where N c = 3 is the color number. 

Due to the spontaneous chiral symmetry breaking, there exist (N 2 — 1) massless Gold- 
stone bosons, where Nf is the light quark flavor number. The S-wave combination Bit 
has the same quantum numbers as the B' Q meson. So the interpolating current ([!]) "sees" 
both B7r and B' . In other words, the contribution due to Bn intermediate states should 
be included explicitly when we write the spectral density at the phenomenological side. 
Otherwise the B' pole contribution will be overestimated leading to an abnormal large 
residue. 

p(s) = \fl,iKs - A +J i) + Ar(a) + • • • , (25) 
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where the first term is the B' Q pole and the p 7r (s) is the Bit intermediate states contribution. 
The excited states and the continuum contribution is denoted by the ellipse. 

We start from the full QCD Lagrangian to derive p^s). It was shown in [pi]] that 



PAP 2 ) = I (0 / ^0(k )6( qo )S(k* - m|)% 2 ) 

(2tt) 3 5(p - k - q)n6(p - k - q ) E I (0|j B , (0)|Btt>| * 
= Ti(l-^)E|(0|j^(0)|i?7r)| 2 

= li(l-^)f(^m/-i) ; (26) 



where the sum is over all the possible Goldstone bosons, = 132MeV, js' = qb is the 
interpolating current of B' in full QCD, f B is defined as, 



< 0\B W q\B(k) >= ifsk^ . (27) 



In the last step of deriving fl26|) the soft pion theorem has been used to calculate the 
matrix element. 

< 0\bq\B(k)n(q) >^ = - < 0\bi l5 q\B(k) >= -if ™ B . (28) 

Moreover the chiral limit m q has been used for all the light quarks and the SU(Nf) 
flavor symmetry has been used to relate the amplitudes for different Goldstone bosons. 
The factor (Nf — -^-) in ( p6|) is the result of summing over these states. 

Letting tub = rrib + A_ i, p 2 = ml + 2rribS and taking the heavy quark limit, i.e., 
nib — ► oo, we have 

J 7T / 



where we have used the relation /_ i = ^/mbfB in the leading order of 1/mQ. Note Ptt(s) is 
a linear function of s while the free parton level spectral density p q (s) in ( pip is quadratic 
in s. The latter is suppressed compared with the former in the lower energy region. So the 
spectral density is significantly disturbed by the presence of light Goldstone bosons. In 
( |29| ) we have introduced the factor #(A + i — s). The reason is that the soft pion theorem 
does not hold any more beyond the region \q\ n < A + i — A_ ( i ~ 350MeV. Moreover it 

was conjectured in [HJ that < 0\bq\B7i > drops when the total energy of Bit becomes 
larger than the mass of B' so that the quark-hadron duality is restored after integrating 
the energy over a larger interval from zero to the continuum threshold. 

Note rrix = 498MeV and m v = 547MeV due to nonzero current quark mass. So in 
realistic case only Bn intermediate states contribute to p n (s) in ( P5|) corresponding to 
Nf = 2. Now we arrive at the new sum rules after making Borel transformation: 

where sq is the continuum threshold. Starting from so we have modeled the phenomeno- 
logical spectral density with the free parton-like one. 
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3.2 New sum rules for g' 

Similarly we can write the double dispersion relation in the leading order of HQET for 

) as 

n(w, J) = - 2 f P( ';f ) dads' + • • • , (31) 

7T J [S — UJ){S' — UJ'J 

where the ellipse denotes the subtraction terms. 
The pole term is 

Pp (s, s') = ^g'U,J_,8{s - k_ tk W - A +) |) . (32) 

and the contribution of the Bn intermediate states in full QCD is 

p^s, s')=J^J <^6{k Q )6(l Q )5{e - ml)5(P) 
(2nf5{p- k - T)n6(p -ko- k) EMq)\jB(0)\B(k)n(l)}(B(kMl)\jB^M\ • ( 33 ) 
Using the soft pion limit and 577/(2) symmetry we find 



3 /b/b'5 1 ' 



PAs, s') = ™* (1 - ^f)8(s - ml,) . (34) 

32 m b + m g s' 



Taking the heavy quark limit (|34[) is reduced to 



Finally we have a new sum rule for 



A i A 



A 



9'L,if + i = 2F n e^{e~^ + ^ f A + J? (s> - A^e-lrds'}- 1 

{-^( Mo )T 2 /i(f ) + MpWT/o(|) + <7i(uo)} • (36) 

3.3 Numerical analysis 

As input we need A_ il / 2 = 0.5 GeV and /_ 1/2 ~ 0.35 GeV 3/2 at the order a s = fl@. The 
numerical results for A +)1 / 2 , f+,1/2 and A +j3 / 2 , /+,3/ 2 in [BJ were obtained by first applying 
the operator ^f^fy to @ and (pTOj ) to extract A +1 / 2 and A + 3 / 2 , which were then used 
to obtain /+,i/ 2 and /+,3/ 2 respectively. Here we use a different procedure which appears 
to be better. This involves with simultaneously varying the parameters A + 1 , / , 1, s etc 
to find the best fitting of the left hand side (L.H.S.) and right hand side (R.H.S.) of the 
sum rules. We work at the region T > 0.4GeV for Eq. ©, where the power correction 
is under control. We allow the continuum threshold to vary from 1.06 GeV to 1.46GeV. 
Numerically we have 

A + ( i = (0.85 ± 0.15)GeV , (37) 

= (0.36 ± 0.10)GeVi . (38) 

With these parameters the left hand side and right hand side agree within five percent in 
the region 0.5 < T < 0.8GeV as can be seen from FIG. 1. Typically at T = 0.4 GeV the 
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sum of the B' Q pole and Bn intermediate states constitutes about 60% of the whole sum 
rule. The continuum starting from sq is about 40%. 

It is important to notice that the Bn intermediate states contribute about 15% to the 
left hand side of (|30|). If we use this fitting method in the numerical analysis of the old 
sum rules (|6]) we reproduce the results in |2j 

A + i = (1.2 ± 0.2)GeV , (39) 

/ +1 = (0.75 ± 0.15)GeVi , (40) 

s = (1.8±0.2)GeV. (41) 

Note in no error is given for / + i . The value and error in (|40|) is the result of our 

reanalysis. We see that both /, i and A + i are significantly reduced after taking into 
account Bn intermediate states. 

We can also apply the fitting method to the analysis of the sum rules for the (1 + , 2 + ) 
doublet. The fitted curve of the L.H.S. and the parton level curve of the R.H.S. of Eq. 
([T0| ) are shown in FIG. 2 in the region T > 0.4GeV with the following most suitable 
parameters A + 3, / + 3, s . 

A +> a = (0.95 ± 0.10)GeV , (42) 



/ + ,| = (0.263 ± 0.06)GeV2 , (43) 
s = (1.3±0.2)GeV. (44) 

With these parameters the fitting is excellent, typically with an accuracy within one 
percent in the large interval 0.5 < T < 1.4GeV as can be seen from FIG. 2. 

The central value of A + 1 is about 100 MeV lower than that of A + 3 , in good agreement 
with the experimental data 



>2 ^'2 



Now we are ready to extract g'. Using the same pion wave functions as in |3 1| we have 



g'f-lf+,1 = (0-36 ± 0.05) GeV 3 , (45) 

where the error refers to the variations with T and Sq. And the central value corresponds 
to T = 0.9GeV and s = 1.26GeV. The variation of the left hand side of ( p6|) with T and 
So is presented in FIG. 2. Finally we get 

g' = 2.8 ± 0.5 . (46) 

The decay width formulas in the leading order of 1 /rag are 

Y{B'^Bn) = ^g^, (47) 

T(B[^B*n) = ^-g' 2 \q\ 7T . (48) 

In order to include the large 1/rriQ correction in the kinematical factors, we use the decay 
width formulas with finite mq instead of (0), (f!8|). 



rW Bx) = ^ % a, , (49) 



7 



1 2 m B .(m%, -m%.f (m 2 B ,+m%.) 2 



'1 "1 



Numerically we have 



T(B' -> fivr) « 250MeV , (51) 
r(Bj -> 5V) » 250MeV , (52) 



with = tub* + 350MeV, mg' = mg + 350MeV. 

4 Summary 

In summary, we have reanalyzed the QCD sum rules for both the (0 + , 1 + ) mass and its 
pionic decay amplitude. The contribution of the Bw intermediate states in the soft pion 
region are investigated in detail. The spectral density is disturbed significantly by the 
presence of Goldstone bosons. After subtracting the contribution of these intermediate 
states, we have obtained new results for the binding energy A + i and decay widths for the 
(0+, 1+) doublet in the leading order of HQET. Numerically A + i = (0.85 ± 0.15)GeV, 
which is about 100 MeV smaller than A + 3 = (0.95 ± 0.10)GeV extracted with the same 
fitting method, in good agreement with the most recent experimental data ||. The 
(0 + , 1 + ) decay width is around 250 MeV, which remains to be larger than the experimental 
result (76 ± 28(stat) ± 15(syst)) MeV in ||. The origin of this discrepancy is not clear at 
present. We want to point out that the same contamination from the Goldstone bosons 
exists for the sum rules for the (1 + , 2 + ) doublet |4], ||. But in this case the Bn intermediate 
states disturb the spectral density only slightly for they are in the D-wave state. 

The errors for our numerical results given in Sec. [O] include only those from the 
variation of the Borel parameter T and the continuum thresold Sq in the window. They 
don't include those from the uncertainty of the condensates and intrinsic errors of the QCD 
sum rule approach. In our analysis we have neglected the contribution of Bit intermediate 
states to the spectral density for s > A + 1 since the soft pion theorem does not hold any 
more and there is not a reliable way to estimate the matrix element < 0\bq\B(k)7c(q) >. 
This is another source of uncertainty. After we submitted the original version of this 
paper, we learned that CLEO collaboration have measured the D' Q mass and width to be 



2461 MeV and 200 ~ 400 MeV respectively |L4|. 
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Figure Captions 



FIG. 1. The varaition of the right and left hand side of Eq. (3C) with Borel 
parameter T is plotted as solid and dotted curves respectively with the fitting 
parameters in (|37l)-(j3§|). 



FIG. 2. The varaition of the right and left hand side of Eq. (|l0|) with T using 
the central values of the fitting parameters in (|42|)-(f44|). 



FIG. 3. The sum rules for g'f_ if, i with sq = 1.36, 1.26, 1.16 GeV respec- 

> 2 ' 2 

tively. 
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